Rapidly variable Fe Ka line in NGC 4051 

J.X. Wangi, Y.Y. Zhou^-^'^, H.G. Xu^, T.G. Wang^'^ 

Received ; accepted 



^Center for Astrophysics, University of Science and Technology of China, Hefei, Anhui, 
230026, P. R. China; jxw@mail.ustc.edu.cn 

^Institute for Space Astrophysics, Department of Physics, Shanghai Jiao Tong University, 
Shanghai, 200030, P. R. China 

^Beijing Astrophysics Center, Beijing, 100080, P. R. China 

■^Laboratory of Cosmic Ray and High Energy Astrophysics, Chinese Academy of Sciences, 
Beijing, 100039, P. R. China 

^The Institute of Physical and Chemical Research (RIKEN), 2-1, Horosawa, Wako, 
Saitama, 351-0198, Japan 



- 2 - 



ABSTRACT 

We present a detailed analysis on the variability of the Fe K emission line in 
NGC 4051 using ASCA data. Through simple Gaussian line fits, we find not 
only obvious Fe K line variability with no significant difference in the X-ray 
continuum flux between two ASCA observations which were separated by ~ 
440 days, but also rapid variability of Fe K line on time scales ~ 10^ s within 
the second observation. During the second observation, the hue is strong (EW 
= 73311?^ eV) and broad {a = 0.96to:i keV) when the source is brightest, 
and become weaker (EW = lOS^gg eV) and narrower (a < 0.09 keV) whilst 
the source is weakest. The equivalent width of Fe K line correlates positively 
with the continuum flux, which shows an opposite trend with another Seyfert 1 
galaxy MCG -6-30-15. 

Subject headings: galaxies: individual (NGC 4051) — galaxies: Seyfert — 
X-rays: galaxies 
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1. Introduction 

It has long been recognized that measurement of the profile of the Fe Ka fluorescence 
line found in many AGNs at ~6.4 keV (Mushotzky 1995, Tanaka et al. 1995, Yaqoob 
et al. 1995, Fabian et al. 1995, Nandra et al. 1997a and references therein) can provide 
an important tracer of matter in the vicinity of the postulated supermassive black hole. 
Doppler and gravitational shifts would imprint characteristic signatures on the line profile 
which map the geometric and dynamical distributions of matter surround the black hole. 
Additional information concerning the geometry of matter in the active nucleus could in 
principle be derived by studying the rapid variability of the line profile, intensity and their 
relationship with the continuum variations. So far, rapid variability in Fe K line has been 
detected only in two Seyfert galaxies. 

Yaqoob el al. (1996) presented the evidence for rapid variability of the Fe K hne 
profile in the narrow-line Seyfert galaxy NGC 7314, which is consistent with a diskline of 
constant equivalent width superposed on a constant flux narrow line (presumably from 
the torus). They found while the X-ray continuum flux varied by a factor of 2 on a time 
scale of hundreds of seconds, the emission in the red wing of the Fe K line below ~6 keV 
responded to those variations on the time scales of less than ~ 3 x 10^ s, and the response 
becomes slower and slower towards the line peak near 6.4 keV. Rapid variability of the Fe 
K line profile was also found in the bright Seyfert 1 galaxy MCG -6-30-15 (Iwasawa et al. 
1996). On time scales of less than ~ 10^ s, the variability pattern is the same as seen in 
NGC 7314. But on long time scales of about several 10^ s, it shows a different behavior. 
When the source is bright, the Fe K hne is weak and dominated by the narrow core, whilst 
during a deep minimum, a huge red tail appears. The intensity of broad Fe K line correlates 
inversely with the continuum flux, so does the equivalent width. 

In this Letter we present an evidence of rapid variability of Fe K line in another nearby 
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(Z = 0.0023) low luminosity narrow-line Seyfert 1 galaxy NGC 4051, which is well known 
for its rapid variability in the X-ray band. Nandra & Pounds (1994) found an equivalent 
width of Fe K line 140 ± 70 eV from the Ginga data by assuming a single narrow line. A 
narrow Fe K emission hne was also detected with an upper limit on the FWHM of ~ 460 
eV and equivalent width EW ~ 170 eV during the first ASCA observation in 1993 (Mihara 
el al. 1994, hereafter M94) for this object. When the target was rc-obscrved by ASCA in 
1994 (Guainazzi el al. 1996, hereafter G96), the Fe K hne was found to be stronger (EW = 
350lJ^g eV) and broad((7 > 0.2keV). 



2. The ASCA data 

NGC 4051 was observed twice by ASCA from 1993 April 25 22:30 to April 26 21:30 
and from 1994 July 6 15:28 to July 9 10:40. In this Letter, wc concentrated on the 
Solid-state Imaging Spectrometer (SIS) data because of the better energy resolution it 
provides (Inoue 1993). The data reductions were performed with the ASCA standard 
software XSELECT according to the following criteria: sateUite not passing through the 
South Atlantic Anomaly, geomagnetic cutoff rigidity greater than 6 GeVc~^, and minimum 
elevation angle above Earth's limb of 10° and 20° for nighttime and daytime observations, 
respectively. Source counts were extracted from a circular area of radius 4' for the SISO and 
SISl. Because of an error in satellite pointing during the first observation, the image fell 
into the dead region between two chips in SISl (M94). Therefore we only used SISO data 
from this observation for spectral analysis. To avoid the complexity caused by the "warm 
absorber" and "soft excess" (M94, G96), we use the 2.5-10.0 keV data for spectral analysis. 
Spectral fits were carried out using the XSPEC, and the background was taken from the 
blank sky data. 
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3. Variability of Iron K Line 

3.1. Comparisons between two observations 

Time average spectra of Fe Ka line profiles for the two observations were given by 
Nandra el al. (1997a) & G96. Due to the poor statistics, the line profile cannot be measured 
exactly, therefore we model the line with a single Gaussian function. The underlying 
continuum is fitted with a single power-law absorbed by the Galactic column density 
(A^^ Ril.3x 10^° cm~^). We did not take into account of the possible Compton refiection 
component here because of its small impact on the ASCA spectra. 

Results of the single Gaussian fits are given in Table 1. For comparison, the results 
of diskline (Fabian el al. 1989, George & Fabian 1991) fits can be found in Nandra el al. 
(1997a) & G96. The variability of Fe K fine is evident (EW changes from 166 eV to 330 eV 
and FWHM from 4500 km/s to 50000 km/s), however, there is no obvious variability in 
the photon index of X-ray continuum. Because of the pointing error of the satellite during 
the first observation, part of the photons fell into the gaps between the CCD chips. The 
continuum fiux of Obsl in Table 1 is only a lower limit. However, since the gaps occupy a 
fairly small fraction of the selected area that includes our target (~ 18%), the continuum 
flux obtained cannot be biased too much. Thus wc conclude that there is no significant 
difference in the continuum fiux between two observations. 



3.2. Rapid FeK Line veiriability in the second observation 

During the two ASCA observations, NGC 4051 shows large amplitude X-ray continuum 
fiux changes on time scales of ~100 s. In order to see if there are short time scale variations 
of Fe K line within a single observation, we choose Obs2 which last longer time (~3 days) 
for analysis in detail. Including a large fiare in the beginning and a deep minimum near the 
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end, the whole observation has been divided into five time intervals (from il to i5) with a 
similar exposure time (~10 ks) as shown in Fig. 1 (diff^erent from G96). Among the five 
intervals, i-2, i-3, i-4 have similar continuum fluxes. 

Results of single Gaussian fittings are shown in Table 1 (We also give the results of 
diskline fits to the five intervals, respectively). The photon index is quite similar between 
i2 to i4 at about F ~1.80. A steeper continuum is suggested for i-1 and a fiatter one for 
i-5. The single Gaussian line dispersions and equivalent widths are also similar between i2 
to i4. A stronger and broader line is detected in i-1, and a weaker narrower one in i-5 (Fig. 
2). In order to show the variability of Fe K line more clearly, we present a contour plot of 
Gaussian line width versus the equivalent width for i-1, i-5 and i-2-|-3-|-4 (a summed dataset 
of i-2, i-3, i-4) obtained from the single Gaussian fits in Figure 3. For i-1, the equivalent 
width of Fe K line is 733 eV, much larger than that for the time average spectrum, and 
the Fe K line is much broader (cr = 0.96 keV). While for i-5, the EW of Fe K line is only 
165 eV, and much narrower (only a narrow fine has been detected, a <0.09 keV). Due to 
the poorer statistics in i-5, there is concerning that a weak broad Fe K component might 
be un-noticeable. To exclude this possibility, we add a broad Gaussian line (Eg — 6.4 keV, 
a =0.46 keV, the same with the time average spectrum) to the fit. The upper limit for the 
EW of this component is 80 eV, far smaller than those in i-1 to i-4. These results clearly 
demonstrate the large variability of Fe K line during the second ASCA observation. G96 
also analyzed the Fe K fine, but failed to detect the variabihty of Fe K line profiles and 
equivalent widths. We would like to point out here that the time sequence selection used in 
G96, different to what we did in this paper, tends to smear out the line variations. 



4. Discussions 
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4.1. The X-ray continuum flux versus the equivalent width of Fe K line 

We detect significant change not only in line equivalent width but also in line profile 
through single Gaussian fitting to the two observations. Prom 1993.4.25 to 1994.7.6, Fe 
K line in NGC 4051 became broader and stronger, in spite of the fact the continuum flux 
remains almost constantly. This shows that the long time scale variability of Fe K line 
can be independent of the continuum variation. Such kind of behavior can be induced 
by the change in structure of the line emission region or the geometry of the hard X-ray 
source in the context of disk line model. The broader and stronger line emission during the 
1994 epoch indicates that the average line emission region is much closer to the putative 
accretion disk. We also detect rapid variability (~ lO^s) of Fe K line within the second 
observation. During the bright flare (i-1), the Fe K line is broad and strong, while during 
the deep minimum (i-5), Fe K line is narrower and weaker. 

Using the ASCA data of 39 AGNs, Nandra el al. (1997b) found a clear decrease in 
the strength of the Fe Ka line with the increase in the luminosity (X-ray Baldwin effect, 
Iwasawa & Taniguchi 1993), which was thought due to the fact that high- luminosity AGNs 
have high accretion rates, causing the accretion disk to become ionized (Matt et al. 1996). 
An inverse correlation between the Fe K line EW and the X-ray continuum flux was also 
found in MCG -6-30-15. The line EWs (sum of the broad and narrow component, calculated 
at the centroid energy 6.4 keV, see Table 2, Iwasawa et al. 1996) versus continuum flux 
is plotted here in Fig. 4b. The case for the NGC 4051 shows that the situation is more 
complicated than this. We find an increase in the equivalent width with no apparent 
change in the continuum flux when the two ASCA observations are considered. Moreover, 
an increase in the line EW with the increase X-ray continuum flux (Fig. 4a), completely 
opposite to the Baldwin effect and the variability in MCG -6-30-15, is found during the 
second ASCA observation. It indicates a complicate mechanism should be considered. 
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4.2. rapid VEiriability on time scale ~ lO^s 

It is well known that the X-ray source in many AGNs is highly variable on short time 
scales ~100 s. As Fe K line is produced via fluorescent process, its strength and profile 
should response to the continuum variability. Obviously, the measurement of the temporal 
response is the best way to map the matter distribution in innermost regions around the 
black hole. In principle, the distance from the X-ray continuum source to the line emission 
region could be determined through observing reverberation effect (Reynolds et al. 1998). 
With the Fe K line profile, we could derive the mass of the black hole in the heart of 
nucleus. The characteristic time scales on which reverberation effects occur is the fight 
crossing time of one gravitational radius, which is ~50 s for a 10^ Mq black hole. But, all 
current available instruments are not able to measure such a fast time response. Even for 
a bright Seyfert galaxy, ASCA can only obtain 1~2 Fe K line photons in 100 s, too few to 
define the fine fiux and profile. However, the Fe K fine variabifity of time scales ~ 10^ s can 
provide constraints on the long time scale (~ 10^ s) variation of the X-ray source or the line 
emission region. So far, such variability has been detected in only three Seyfert galaxies, 
NGC 7314, MCG -6-30-15 and NGC 4051. For NGC 7314, a simple disk-plus-torus model 
could explain its rapid variability fairly well. But for MCG -6-30-15 and NGC 4051, the 
situation is more complicated. 

According to the disk model, the fiuorescent iron fines is produced from the hard X-ray 
irradiation of the disk composed of cold gas around a black hole. When the X-ray source is 
bright, iron will start to be ionized, and in the intermediate ionization states, the resonance 
scattering can cause a reduction in the line fiux (Matt et al. 1996). It is a possible reason 
causing the anti-correlation of the EW of broad components to the X-ray continuum flux in 
MCG -6-30-15. We consider that the ionized disk is also a possible reason for the positive 
correlation in NGC 4051 because a factor ~2 larger EW for FeXXV than cold iron could 
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occur (Matt et al. 1996). Though not well constrained, the diskhne energies for il to i5 in 
Table 1. show an trend to correlate with the continuum flux, which gives some support to 
this opinion. 

There is also another possible explanation to the the rapid variability of Fe K line. 
Fabian (1997) suggested that the X-ray continuum might be generated by magnetic flares 
above the accretion disk. At any given time there are only a few flares otherwise the 
rapid variability will be average out. When the dominated flares move around on the disk, 
changes in the Fe K line profile would be seen. Such model could explain the odd variability 
of Fe K line in MCG -6-30-15 well. We think it is also appropriate for the rapid variability 
in NGC 4051 discovered in this paper. For the sequence i-5, the Fe K fine profile is similar 
to the bright fiare in MCG -6-30-15. The small FWHM might due to a dominate fiare 
located far away from the central black hole, or a succession of fiares on the approaching 
side of the disk. The broad line during i-1 might due to a bright dominate fiare very close 
(~6i?p) to the black hole, and the large EW might due to a overabundance of iron. 

This work is supported by Chinese National Natural Science Foundation, PanDeng 
Project and Foundation of Ministry of Education. The authors would also hke to thank 
S.A. Huang for great help with ASCA data reduction. 
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Observing Time: 1994.7.6 Bin Time: 128.0 s 




5x10^ 10^ 1.5x10^ 

Time (s) 



Fig. 1. — The 0.4-10.0 keV light curve from the SISO for the second ASCA observation, 
binned at 128 s. Time-intervals used in the section 3.2 are also indicated here. 
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Fig. 2. — Ratios of the data to an absorbed power-law model fit to the data in the 2.5 to 
4.5 keV and 7.0 to 10.0 keV bands to illustrate the Fe K line profile for the SIS in the time 
intervals i-l,i-2-|-3+4,i-5, which are indicated in Fig. 1. 




Fig. 3. — 68 and 90 per cent contours in the line width and the hne EW for two degrees 
freedom from the single Gaussian fit to the spectra in the time-intervals i-1, 1-2+3+4, i-5. 
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Fig. 4.— Plot of Fe K line EWs versus continuum flux for: a) NGC 4051 and b)MCG 
-6-30-15. 
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Table 1: Single Gaussian fits to the two ASCA observations and the five datasets of Obs2 





pa 


Eg (keV) 


a (keV) 


EW (eV)^ 


x'/dof 


flux'' (keV) 


Ri/Rg 


Obsl 








162t?i 


148/127 


> 2.15 




Obs2 


1 oq+0.05 
1.00_o.06 




0.46lg:}! 


3301^^ 


308/335 


2.42 




i-1 




6.451^:^ 


0.96lg:i 


70 0+206 
' •^•J-219 


184/190 


... 6.74j:gig 


6.31^:1 


i-2 


1 71+0.14 

' -L-o.is 


6.0211 


U-"±"±-0.18 


2961}^^ 


174/172 


... 6.491°;}} 


13.1^}- 


i-3 


1 oo+O.lO 
i.»^_0.10 


fi 47+0-19 

-0.13 


0.211?^ 


269li 


228/227 


... 6.49lgi} 




i-4 






0.46l°i^ 


453l}«^ 


119/138 


... QMtlil 




i-5 


1 CI +0.15 
J^-<3i-0.15 






16511^ 


152/149 


... 6.401H^ 


9991^87 



"Best fitting value when data in 4.5-7.0 keV region are excluded 

X 10-^^erg.cm-2.s-^ (2-10 keV) 
'^calculated at the centroid energy 6.4 keV 

'^Results of diskhne models fits with disk emissivity index q fixed at -3.0, the outer disk radii 
Ro at 1000i?g {Rg = GM/c^) and the inclination of the disk at 25° (G96) 



